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ZNTRCDUCTICN 


In  the  prediction  of  transients  involving  hydrodynamic  or  aerodynamic  phe¬ 
nomena*  it  has  been  customary  to  calculate  pressure  variation  and  fluid  resistance 
neglecting  possible  effects  ef  unsteadiness  on  the  mechanics  of  the  fluid  motion. 

There  are  several  areas  for  which  knowledge  of  the  effects  of  unsteadiness  would  be 
useful.  Included*  for  example*  are  transient  resistance  and  stability  of  accelera¬ 
ting  missiles  and  other  immersed  objects*  flow  meter  coefficients  with  pulsating 
flows*  and  transient  hydrodynamic  performance  of  pumps*  compressors  and  turbines. 

The  latter  is  a  case  where  under  many  circumstances  steady  state  performance  has 
been  used  with  good  results  in  predicting  transient  pressures  and  machine  accelera¬ 
tions  and  decelerations.  Yet*  recently*  in  connection  with  pumps  of  special  design* 
discrepancies  between  measured  and  calculated  transients  have  indicated  what  is 
probably  an  effect  of  unsteadiness  on  the  basic  fluid  motion.  In  all  cases  a  basic 
question  is  the  effect  of  unsteadiness  on  fluid  shear  and  turbulence  generation*  and  4 
the  resulting  effects  on  the  inertial  and  frictional  components  contributing  to  the  a 
instantaneous  total  potential  drop.  1 

This  paper  summarizes  the  results  of  investigations  in  the  M.  I.  T.  Unsteady 
Flow  Water  Tunnel  (Ref.  1*  2*  3*  4)  of  accelerated  and  decelerated  flow  through  uni¬ 
form  conduits  and  orifices  in  conduits.  In  the  uniform  conduit  shear  and  turbulence 
is  generated  through  boundary  layer  friction  and  is  essentially  unifozm  along  the 
duct.  The  orifices  cause  separation  and  jet  formation  with  accompanying  high  shear  3 
and  turbulence  which  varies  along  the  duct  as  the  jet  diffuses  and  the  turbulence 
is  dissipated.  I  21st  I 


UNSTEADY  FLOW  EQUATIONS 


MM 


As  a  basis  for  analyzing  experimental  results*  the  following  momentum 
analysis  was  made  in  which  the  effects  of  all  the  variables  including  turbulence 
are  considered.  Thus  while  the  equations  are  reduced  ultimately  to  essentially 
a  one-dimensional  form*  more  insight  is  given  to  the  significance  of  each  term 
than  is  readily  apparent  from  an  ordinary  one-dimensional  analysis. 

Consider  a  constant  diameter  conduit  which  may  be  unobstructed  or  may 
contain  constrictions  such  as  orifices  or  venturi  sections.  Let  the  dotted  boundary 
shown  in  Figure  1  define  a  control  volume  for  a  general  case. 
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2. 


Applying  the  momentum  principle  to  unsteady  turbulent  liquid  flow  through  such  a 
system#  we  can  write  the  instantaneous  balance 

(External  Forces)  *  (Net  flux  of  momentum  from  control  volume  +  rate 

of  change  of  momentum  within  the  volume) 

or#  for  the  x-direction 

rA  l  rA  2  rA2  r*  l  f * 

J  PxdA  -  J  P2dA  -  F  =  /ru22dA  -  (l) 

where  P  =  x-component  of  local  pressure  intensity 

F  =  lumped  boundary  resistance  due  to  wall  shear  and 
any  constrictions 

U  =  local  velocity 

¥  =  liquid  volume  within  control  surface 

Analogous  to  "steady"  turbulent  flow  we  will  assume  that  at  each  instant 
during  a  transient  the  velocity  field  can  be  described  as  a  "mean  flow"  plus  a 
"fluctuating  flow"  and  introduce  the  usual  notation 

U  =  u  +  u«  (2) 


Here  then  u  is  a  kind  of  mean  value  like  the  temporal  mean  value  of  statistically 
steady  turbulent  motion.  Using  Eq,  (2)  gives 
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(3) 


At  this  stage  it  is  customary  in  analyzing  turbulent  flows  to  take  temporal 
mean  values  of  all  terms  in  the  equations  in  order  to  obtain  a  relation  governing  the 
main  flow.  This  ignores  the  contribution  of  the  terms  that  are  linear  in  the  tur¬ 
bulent  components#  because  they  drop  out  in  the  averaging  process.  For  very  rapid 
changes  of  the  main  flow  this  may  not  be  justified.  We  will  leave  the  equation  in 
terms  of  the  instantaneous  velocity  fluctuations  so  that  the  relative  importance  of 
the  instantaneous  turbulence  can  be  evaluated.  Using  the  definitions) 


IfdA  _ 
p_  - 
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average  pressure  intensity  over  the  cross  section 
instantaneous  cross-sectional  mean  velocity  in  conduit 

unsteady  flow  coefficient  of  total  drop  in  potential 


(4) 


3. 
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unsteady  flow  coefficient  of  boundary  resistance 


distribution  factor  for  mean  velocity 

r , 

&  . ,  HU  dA  =  distribution  factor  for  turbulent  fluctuations 
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and  the  approximation  that  ^  is  independent  of  time*  Eq.  (3)  is  expressed  dimen* 
sionlessly  as  uo 


K  =  K  +  2 
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=  acceleration  of  the  mean  flow  in  the  conduit 
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The  last  term  in  Eq.  (8)  is  the  dimensionless  force  of  inertia  of  the 
turbulent  fluid  to  local  accelerations.  The  second  right-hand  texm  gives  the  ef¬ 
fect  of  non-uni foxmity  of  the  mean  flow  and  turbulence  intensity  between  sections 
1  and  2.  It  is  the  dimensionless  x  -  component  of  the  flux  of  momentum  of  the 
absolute  motion.  While  the  effect  of  turbulence  appears  explicitly  only  in  the 
second  and  third  right-hand  terms*  it  also  appears  indirectly  in  the  resistance 
coefficient  K.  The  velocity  and  turbulence  distributions  within  the  liquid  volume 
are  interdependent  with  the  boundary  shears  and  pressures  and  hence  the  boundary 
resistance. 


Also  to  the  extent  that  the  establishment  of  each  instantaneous  velocity 
and  turbulence  distribution  requires  some  absolute  time  interval  for  the  adjustment 
from  a  previous  condition*  the  relative  magnitudes  of  the  several  terms  in  this  force- 
momentum  balance  may  depend  on  the  rate  of  change  of  acceleration  as  well  as  the  mag¬ 
nitude  of  the  acceleration. 
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While  Eq.  (8)  is  useful  for  qualitative  indication  of  how  turbulence* 
flow  uniformity  and  acceleration  affect  the  total  potential  drop  and  boundary  re* 
si  stance,  it  cannot  be  used  for  quantitative  comparisons  because  the  retired 
instantaneous  velocity  and  turbulence  distributions  cannot  be  calculated  or  measured. 
Because  of  this  and  because  it  is  desirable  to  make  comparisons  with  steady-state 
conditions  corresponding  to  given  instantaneous  rate  of  discharge  we  introduce  the 
following  simplifications! 

K  (u) 

=  K  +  K. 


K  =  unsteady  flow  coefficient  of  boundary  resistance  and 
u  momentum  flux  of  absolute  local  velocity 

K  =  one-dimensional  steady  state  "resistance"  coefficient 

5 

Kt  =  correcting  coefficient  to  measure  the  additional  trans¬ 
ient  effects  on  boundary  resistance  and  momentum  flux  of 
absolute  local  velocity 

Eq.  (8)  becomes 

K  =  K  +  K  +  c.  ^  (12) 

a  s  t  1  u  2A 

o 

For  steady  flow  this  reduces  to  the  relation  given  by  the  conventional  one-dimensional 
energy  equation  where  the  resistance  coefficient  is  taken  as  a  measure  of  the  energy 
dissipation. 

Eq.  ( 12)  can  be  simplified  further  with  the  aid  of  an  analogy  to  Schonfeld's 
analysis  for  smooth  round  tubes  (Ref.  5).  Schonfeld  presents  the  following  solution 
for  the  special  case  of  slowly  varied  motion  in  which  the  resistance  dominates  (as 
opposed  to  quickly  varied  motion  where  the  inertia  dominates). 


Lets 


where 


where 


Pi 
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dQ 

dt 


Q  =  rate  of  discharge 
R^  =  hydraulic  radius 
C1  =  steady  flow  Chezy  coefficient" 

N  =  -7^  Ci.0  - - ^ 

*  (c*  +  u.or 

By  substituting  C*  (fg  =  steady  flew  friction  factor), 

R^  =  ^  and  Q  =  AUq,  Eq.  (14)  can  be  rewritten  for  the  tunnel  test  section  thus 
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5. 


We  note  that  this  can  be  put  in  the  form  of  Eq.  (2)  if  we  make  the  substitutions 
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Cj  =  1.00 
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(17) 


As  these  substitutions  indicate*  Schonfe Id’s  analysis  considers  the  turbulence  to 
be  fully  developed  at  every  instant  and  uniform  conditions  to  exist  along  the  length 
of  the  conduit.  Thus  the  value  c^  =  1.00  implies  that  the  mean  value  of  the  turbu¬ 
lent  fluctuations  over  the  liquid  volume  of  the  uniform  tube  is  zero  at  all  times 
so  that  the  inertial  coefficient  is  merely 

'  <■-  dV 
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Furthermore*  if  similar  flow  conditions  exist  at  all  sections  along  the  conduit* 
the  term  (<tf"2  -oCj)  will  be  zero  and  the  term  (l2  -  Ij)  will  approach  zero  for 

sufficiently  random  values  of  u*  at  each  instant  over  each  annular  increment  of 
flow  area  at  the  two  sections.  Thus  as  given  by  Eq.  ( 17)  implies  no  influence 

of  non-uniformity  in  velocity  or  turbulence. 


Returning  now  to  the  more  general  case  which  may  include  tubes  with 
constrictions*  let  us  write  in  analogy  to  Eq,  (17) 


Kt  ~  c2 


2a  ■» 
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(18) 


In  Eq.  (18)  c2  is  a  measure  of  the  deviation^  due  to  unsteadiness  of  the  boundary 
resistance  and  flux  of  momentum.  Using  Eq.  (18)*  Eq.  (12)  reduces  tot 
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We  can  also  write  this  as 


c  =  cx  +  c2 


K 

-2 

K 


=  1  ♦  c 


2n9 

— sr 

KFA 
s  o 


(20) 


or*  noting  from  Eqs.  (ll)  and  (12)  that 

K  =  K  +  K  =  K  -  c 
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6. 


Equations  (19)*  (20)  and  (21)  are  exact  statements  of  the  balance  of 
forces  specified  by  the  momentum  principle.  Moreover*  to  the  extent  that  the 
quantity  2^(o<2  “  ^j)  +  "  IiHin  Eq#  approaches  zero,  Ky  becomes  equal 

to  K  and  the  equations  will  be  useful  for  comparing  the  steady  and  unsteady  boundary 
resistance  at  given  instantaneous  discharge  rates.  Here  we  note*  that  as  for  a  clear 
conduit*  flow  through  a  constricted  conduit  having  a  test  length  L  much  longer  than 
the  disturbed  flow  zone  caused  by  the  constriction*  should  exhibit  approximately  the 
same  velocity  and  turbulence  distributions  at  sections  1  and  2.  Hence*  only  a  small 
error  is  introduced  by  assuming  K  =  K,  In  this  case*  it  will  be  noted  from  Eq.  (21) 
that  for  accelerated  flow  where  auis  positive,  positive  values  of  c2  will  indicate 
that  more  boundary  resistance  is  developed  than  for  steady  flow  and  vice  versa. 

For  decelerated  flow  where  a  is  negative,  positive  c2  will  indicate  that  less 

boundary  resistance  is  developed  than  for  steady  flow  and  vice  versa. 

Note  that  in  c^,  as  defined  by  Eq,  (10),  the  magnitude  of  the  first  term 

will  be  unity  or  larger  unless  there  are  flow  reversals  giving  negative  u  values. 
Again  as  for  the  clear  conduit,  the  mean  value  of  the  turbulent  fluctuations  over 
the  liquid  volume  of  the  constricted  conduit  is  expected  to  be  very  small  if  not 
zero*  making  the  second  term  of  c.  negligible.  Let  us  introduoe  then  the  approxi- 
mation  for  c^  rL 

C1^L  I  u“  dx  (22) 

J  O 

o 

where  V  =  average  velocity  over  the  cross  section  of  the  main  stream  (jet)  at  any  x. 

This  definition  peimits  evaluation  of  Cj  from  a  flow  net  of  the  jet  profile  through 
the  constriction.  By  this  definition  also  Cj  will  tend  toward  a  constant  value  if 

the  dimensionless  velocity  along  the  main  stream  (or  jet)  remains  independent  of 
the  unsteadiness.  An  experimental  determination  of  the  coefficient  c  and  a  calcu¬ 
lated  value  of  Cj  will  permit  evaluation  of  c2  and  .  With  these  simpli¬ 

fications*  the  resulting  c2  will  absorb  the  difference  between  the  true  inertial 

coefficient  and  the  value  calculated  by  Eq.  (22),  as  well  as  the  durations*  due 
to  unsteadiness*  of  the  boundary  resistance  and  flux  of  momentum. 

The  parametric  form  of  Eqs.  (20)  and  ( 21 )  makes  it  convenient  to  evalu¬ 
ate  the  comparison  between  steady  and  unsteady  behavior  from  a  simple  series  of 
measurements  of  total  potential  drop  along  the  conduit  versus  instantaneous  flow 
rate.  From  each  such  basic  experiment,  the  comparison  can  be  obtained  for  each  of 
a  range  of  values  of  the  ratio  i £  =  £L.  « 

U  2A  U  2 
o  o 

•vU  gL 

Finally  note  that  the  term  a~r~  =  is  a  parameter  proportional  to  the 

U  ^A  U  * 
o  o 

ratio  of  local  to  convective  acceleration.  If  all  the  effects  of  unsteadiness  are 
fundamentally  dependent  on  the  acceleration,  the  coefficients  in  Eqs.  (20  and  (21) 
will  be  constants,  otherwise  not. 
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EXPERIMENTAL  APPARATUS 

The  apparatus  used  for  these  experiments  is  s  non-return  unsteady  flow 
water  tunnel  (Ref.  6).  As  shown  by  the  schematic  section  in  Figure  2,  the  tunnel 
consists  of  two  cylindrical  tanks  mounted  one  above  the  other  and  connected  by  a  1- 
inch  diameter  smooth  brass  conduit  99  diameters  in  length.  This  conduit  constitutes 
the  test  section  in  which  constrictions  such  as  orifices  or  venturi  sections  can  be 
placed.  Water  is  caused  to  flow  from  one  tank  to  the  other  under  pneumatic  control. 
Compressed  air  is  admitted  to  the  spaces  above  the  water  surfaces  in  the  two  tanks 
to  provide  a  driving  force  for  a  desired  flow  rate  and  acceleration  or  deceleration. 

To  obtain  the  desired  ranges  of  acceleration  and  pressure  in  the  working  section,  com¬ 
pressed  air  must  be  admitted  to  or  released  from  either  tank  according  to  some  time 
schedule.  To  prevent  cavitation  and  the  introduction  of  air  into  the  piezometric 
system,  the  test  section  is  maintained  at  positive  pressure  by  throttling  the  exhaust 
from  the  bottom  tank.  The  square  edged  orifices  employed  were  dimensioned  according 
to  ASME  standards  and  constructed  from  0.102'1  thick  sheet  brass.  The  conduit  and 
orifice  combinations  used  and  the  principal  dimensional  data  and  location  of  piezo¬ 
meter  taps  are  given  in  Table  I  and  Figure  3. 

Table  I 


Test  Combinations  and  Principal  Dimensions 
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Figure  3 

In  all  experiments  the  beginning  of  the  test  length  over  which  potential 
drops  were  measured  was  located  well  beyond  the  calculated  transition  distance 


T 


1 


.. ;  \ 

necessary  to  establish  a  fully  developed  turbulent  .velocity  profile  in  the  1-inch 

tube.  In  most  tests  the  test  length  began  at  approximately  33  diameters  from  the 

entrance  nozzle,  the  orifice  being  located  at  nx'*  =  37.5  diameters.  (See  Figure  2) 

In  some  cases  a  28  diameter  approach  was  used  ("x"  =  32.5  diameters).  The  actual 

test  length  L  was  27  diameters  for  the  unobstructed  tube  and  12  diameters  for  the 

cases  of  orifice  constricted 'tube.  For  the  latter,  the  piezometer  taps  were  located 

U.5  diameters  upstream  and  7.5  diameters  downstream  of  the  orifice  plate.  The  7.5 

diameter  downstream  distance  was  chosen  to  include  the  expected  zone  of  influence  of 

the  orifice  on  the  local  flow  conditions  (Ref.  7).  Some  measurements  were  made  also 

with  the  downstream  piezometer  taps  22.5  diameters  from  the  orifice  as  checks  in 

case  the  orifice  disturbance  persisted  for  greater  distances  for  unsteady  flow  than 

for  steady.  The  results  of  these  checks  did  not  alter  the  conclusions  drawn  from 

the  measurements  over  the  shorter  distance. 

The  nozzle  at  the  inlet  to  the  working  section  is  used  for  flow  measure¬ 
ments.  The  instantaneous  pressure  drops  recorded  during  unsteady  flows  are  corrected 
to  account  for  the  inertia  force  due  to  the  local  acceleration  of  the  fluid  through 
the  nozzle.  The  correction  calculated  assuming  potential  flow  0.2  dUQ/dt. 

The  several  differential  pressures  were  measured  with  diaphragm  type 
pressure  cells  in  which  the  diaphragm  deflection  actuates  a  differential  transformer. 
Each  pressure  gage  signal  is  sent  through  a  separate  amplifying  and  detecting  unit 
and  is  then  recorded  versus  time.  Two  amplifying-recording  systems  were  used 
during  the  experiments;  in  one,  a  photographic  record  was  obtained  using  a  Hathaway 
Type  S8-C  oscillograph;  in  the  other,  the  record  was  traced  with  a  temperature- 
stylus  on  Sanborn  "Permapaper"  using  a  four-channel  Sanborn  Recorder,  Model  150. 
Measured  natural  frequencies  of  the  pressure  cells  and  oscillograph  recording  system 
connected  as  for  testing  including  water-filled  lead  lines,  exceeded  165  cycles  per 
second.  With  the  Sanborn  recorder  the  limit  of  accurate  response  is  about  100  cycles 
per  second.  Instantaneous  differential  pressures  were  evaluated  using  static  cali¬ 
brations  of  each  gage  before  and  after  test  runs. 

EXPERIMENTAL  PROGRAM 

•  The  range  of  variables  for  the  experiments  reported  here  are  given  ir 
Table  II.  As  the  table  indicates,  a  range  of  velocities  and  accelerations  or  de¬ 
celerations  were  included.  In  addition,  the  several  tests  included  different  rates 
of  change  of  acceleration. 

Acceleration  tests  were  made  from  zero  velocity  or  wer^  preceded  by  an 
initial  period  of  steady  flow.  All  deceleration  tests  were  preceded  by  an  initial 
period  of  steady  flow.  At  the  start  of  each  unsteady  period  was  an  initial  impulse 

phase  during  which  the  acceleration  (cr  deceleration)  changed  rapidly  (j|^»  0). 

Following  this  was  an  '’established'*  phase,  distinguished  ty  either  constant  or  more 
slowly  changing  acceleration.  In  general,  the  acceleration  or  deceleration  varied 
continuously  although  in  the  case  of  the  uniform  diameter  conduit,  several  runs, 
each  with  essentially  constant  deceleration,  were  obtained.  The  period  of  the 
initial  impulse  phase  varied  and  it  was  not  always  possible  to  obtain  reliable 
recordings  of  the  instantaneous  pressures.  Typical  test  results  appear  in  Figures 
8a  and  8b.  These  diagrams  show  instantaneous  conduit  velocity  UQ,  acceleration 

dUQ/dt  and  total  potential  head  drop  Ha  as  calculated  from  oscillograph  recordings. 
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Table  II 


Range  of  Experimental  Investigation 


Equiva¬ 

lent 

area 

ratio 

Maximum 
conduit 
velocity 
at  steady 
state 
discharge 

fps 

1 

■ 

Maximum 

conduit 

Reynolds 

Number 

Maximum 

throat 

Reynolds 

Number 

Conduit 

accele¬ 

ration 

range 

fps* 

max.  ave. 

i 

Conduit 

decele¬ 

ration 

range 

fps2 

max.  ave. 

Orifices: 

0.7 

38.7 

320,000 

580,000 

80 

Uo 

50  25 

0.5 

29.75 

21*6,000 

720,000 

60 

30 

50  25 

0.3 

18.00 

150,000 

770,000 

30 

15 

30  15 

Smooth  Conduit: 

1.0 

9  to  60* 

500,000 

— 

80 

Uo 

|  36  to  18** 

1 

300,000 

. 

7,  11, 

16,  20 

*  Runs  39,  UO,  by  Deemer  for  accelerated  flow. 
**  Runs  UJ,  by  Jordaan  for  decelerated  flow. 


H 

From  such  data  values  of  K  =  — 

a  ^ 


versus 


aL 


were  obtained  for  suc- 


o  '  °  o 

cessive  time  intervals  throughout  the  test.  Any  one  run  gives  a  wide  range  of 

values  of  — , 0  .  From  steady  flow  data,  values  of  K  versus  Reynolds  number  were 

U  3 

o. 

computed.  These  were  used  to  relate  the  successive  instantaneous  conditions  of  un¬ 
steady  flow  to  steady  state  conditions  for  the  same  Reynolds  number.  With  computed 
values  of  the  inertial  coefficient  c-^  the  combined  boundary  resistance  and  momentun 

flux  during  unsteady  flow  was  evaluated  using  Eq.  (21). 


As  previously  mentioned,  the  inertial  coefficients  c^  were  calculated  for 

the  orifices  by  numerical  integration  of  Eq.  (22)  from  a  flow  net  of  the  jet  profile. 
This  profile  is  known  only  approximately  and  in  addition  is  assumed  to  be  essentially 
the  same  over  the  range  of  velocity  and  acceleration  of  the  tests.  Therefore,  values 
of  Ku  can  be  determined  only  within  some  range.  Table  III  gives  the  computed  magni¬ 
tudes  of  c^  together  with  extreme  limits  of  possible  deviations. 
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Table  HI 

Inertia  Coefficients  for  Unsteady  Flow  through  Orifices 


Orifice 

area 

ratio 

c,  Based 
x  on  U  dia. 
jet 

diffusion 

length 

c.  Possible 

1  range  of 
variation 

c,  Assumed 

for  i 

calculation 

purposes 

0.7 

l.lli 

1.00  to  1.5c 

. 

1.15 

0  .5 

1.27 

1.00  to  1.80 

1.30 

0.3 

1.51 

1.00  to  2.00 

1.50 

RESULTS  AND  CONCLUSIONS 


Steady  Flow 


The  experimentally  determined  steady  flow  resistance  coefficients  K  are 

s 

presented  in  Table  IV.  The  coefficients  for  orifices  are  in  fair  agreement  with 
values  calculated  from  the  sudden  expansion  formula.  Some  dependence  on  Reynolds 
number  is  indicated  for  the  0.7  area  ratio  unit.  Otherwise  the  coefficients  are 
essentially  constant  over  the  velocity  range  covered  by  the  experiments.  The  pipe 
friction  coefficients  are  given  by  the  relation 

-j=r  -  ©-•&>  '  °-8 

Ts 

U  D 

where  =  Conduit  Reynolds  number  =  —— 

f  =  Steady  state  pipe  friction  factor 
s 

Substituting  Kg  =  fg  g  with  g  =  12  gives  the  equation  in  Table  IV  for  a  foot  length 
of  test  conduit. 


Table  IV 

Steady  State  Discharge  Coefficients 


j  Area  ratio  .  i 

.  ...  .  i  .  _ 

.  Loss  ccef 

s 

Remarks 

Orifices  j  0.7 

0;93 

Low  Velocities 

1 

i 

0.98 

High  Velocities 

! 

3.81 

j  0.3 

! 

17.00 

_ 

. 

Smooth  Conduit  j  1.0 

! 

1A/K^  =  o.59  iog10  (ff-jicy  -  o.5U 

J 
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Unsteady  Flov 

The  comparison  of  the  unsteady  and  steady  behaviors  are  given  in  Figures 

JC  .  • 

U-7  in  the  form  of  diagrams  of  ,  versus  —  ■»  .  Each  point  plotted  on  these 

s  K.Uo 

..so 

diagrams  was  evaluated  using  Kb  from  Table  IV  and  the  inertial  coefficient  from 
Table  III.  8 

First  it  is  noted  that  the  plotted  points  are  spread  over. a  considerable 
area  in  each  diagram.  This  is  due  in  part  to  the  fact  that  this  form  of  represen¬ 
tation  is  sensitive  to  small  differences.  Therefore,  errors  are  exaggerated.  In 
addition,  s  the  data  included  a  range  of  velocities,  accelerations  and  rates  of 
change  c.  acceleration  (or  decelerations)  the  plotted  spread  is  an  indication  that 
the  actual  unsteady  velocity  and  turbulence  distribution  and  resulting  boundary 
shear  and  boundary  pressures  are  in  some  way  dependent  on  these  factors. 

Nevertheless,  it  is  seen  that  in  each  diagram  the  data  falls  essentially 
in  two  opposite  quadrants,  indicating  definite,  even  though  qualitative,  trends  in 
the  relative  magnitudes  of  K  and  K  .  Merely  to  emphasize  these  trends  the  data  is 
represented  by  single  straight  lines  with  positive  or  negative  slopes.  These  slopes 
are  measures  of  Cg  in  Eq.  (21)  and  values  of  c«  are  indicated.  However,  it  is  em¬ 
phasized  that  only  the  sign  of  the  slope,  and  of  c2,  is  significant,  not  the  magni¬ 
tude.  In  drawing  these  lines,  emphasis  was  given  to  the  "established  phase”  portion 
of  the  test  run,  where  the  rate  of  change  of  acceleration  is  not  large. 

For  accelerated  flow  through  the  three  orifices,  the  unsteady  coefficient 
K„  is  less  than  the  steady  K  at  the  same  instantaneous  velocity.  Assuming  as  pre- 

viously  mentioned  that  the  net  flux  of  momentum  is  zero  for  the  volume  between  the 
measuring  stations,  the  indication  is  that  the  boundary  resistance  during  accelera¬ 
tion  is  less  than  for  the  equivalent  steady  motion.  For  decelerated  flows  through 

these  orifices,  K  >  K  and  the  boundary  resistance  exceeds  that  for  steady  flow, 
u  s 

It  is  recognized  that  the  results  are  qualitative,  and,  because  of  the 
uncertainty  in  the  value  of  the  calculated  inertial  coefficient,  this  would  be  the 
case  even  in  the  absence  of  the  factors  which  were  just  mentioned  as  contributing 
to  the  spread  of  observed  data. 

As  previously  described,  the  values  of  c^  used  to  obtain  Cg  and  Ku  were 

calculated  by  numerical  integration  of  Eq.  (22)  from  a  flow  net  of  the  jet  profile. 
This  profile  is  known  only  approximately  and  in  addition  is  assumed  to  be  essentially 
the  same  over  the  range  of  velocities  and  accelerations  of  the  tests.  On  the  other 
hand,  it  should  be  emphasized  that  the  errors  probable  or  possible  in  the  experimental 
measurements,  or  in  determining  the  inertial  head  drop  term,  would  not  alter  the 
stated  conclusions. 

In  the  case  of  the  uniform  diameter  conduit,  the  magnitude  of  the  boundary 

resistance  during  accelerated  flow  is  very  nearly  the  same  as  for  the  equivalent 

steady  motion.  Nevertheless,  in  Figure  7  there  is  a  definite  indication  that  K A 

u  s 

is  greater  than  unity.  Thus  the  case  of  resistance  due  to  boundary  layer  shear 
stresses  is  affected  differently  by  unsteadiness  than  resistance  associated  with  the 
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turbulence  generation  and  diffusion  accompanying  .  separation  and  jet  formation.  For 

these  data,  Schonfeld* s  theory  is  used  as  a  guide,  and  a  straight  line  having  a  small 

positive  slope  is  drawn  through  the  plotted  points.  The  relation  in  Table  IV  shows 

K  to  be  a  function  of  Reynolds  number  while  the  relations  derived  from  Schonfeld*  s 
s 

theory  (Eqs.  15-17)  predicts  Cg  and  Ku  also  to  be  dependent  on  R.  The  variation  in 
Cg  is  small,  however.  Using  Eqs.  (17) .and  (18)  over  the  range  of  Reynolds  number 
investigated,  Cg  varies  only  between  0.010  and  0.015.  Hence  a  single  straight  line 
with  a  slope  indicating. a  constant  Cg  =  0.010  was  arbitrarily  chosen  to  qualitatively 
represent  the  test  data. 

For  decelerated  flow  the  boundary  resistance  of  the  uniform  tube  is  less 
than  for  steady  flow.  In  this  case,  however,  there  is  a  clear  indication  of  effects 
not  predicted  by  Schonfeld* s  results.  As  shown,  these  data  can  be  represented  by  a 
family  of  lines,  essentially  parallel,  one  for  each  deceleration.  At  any  particular 
velocity,  the  proportion  of  boundary  resistance  to  overall  potential  drop  is  differ¬ 
ent,  decreasing  with  increasing  deceleration.  All  of  these  runs  were  started  from 
the  same  steady  state  velocity,  but  included  different  initial  impulse  periods. 

From  the  parallel  displacement  of  the  lines  for  different  decelerations,  it  appears 
that  the  flow  conditions  of  the  subsequent  established  phase  depend  on  the  previous 
flow  history. 

These  observations  for  the  uniform  tube  are  consistent  with  the  view  that 
under  acceleration  the  central  portion  of  the  stream  moves  somewhat  bodily  while  the 
velocity  profile  steepens,  giving  higher  shear.  For  deceleration,  the  reverse  seems 
to  hold.  In  either  event,  it  appears  that ■ unsteadiness  does  not  result  in  marked 
changes  from  equivalent  steady  state  flows. 

In  the  case  of  the  orifices,  however,  it  appears  that  the  imposition  of  a 
transient  results  in  flows  having  quite  different  velocity  and  turbulence  character¬ 
istics.  This  was  indicated  not  only  by  the  relative  magnitude  of  and  Kg,  but 

also  by  what  was  first  thought  to  be  an  anomalous  experimental  result.  For  decele¬ 
rated  flow  through  the  smaller  orifice's,  it  was  observed  that  as  the  unsteady  run 
proceeded  the  magnitude  of  the  potential  drop  changed  from  less  than  the  equivalent 
steady  state  drop  (as. required  to  establish  the  deceleration)  to  more;  i.e.,  K  /K 

cl  S 

became  greater  than  1.0  as  the  test  run  proceeded.  This  observation  was  repeated  on 

many  runs  and  cannot  be  attributed  to  measurement  errors.  For  acceleration  through 

the  0.3  orifice,  there  was  some  indication  that  a  corresponding  change  to  K  ^  1.0 

3  s 

occurred  late  in  the  run,  however, experimental  errors  could  conceivably  account  for 
the  shift  in  this  case.  Such  results  could  only  mean  that  as  the  unsteady  flow  pro¬ 
ceeded  the  internal  structure  of  the  velocity  and  turbulence  distribution  changed 
to  the  point  that  it  was  no  longer  comparable  to  ary  steacfy  state  flow  condition. 

Such  effects  as  mentioned  in  the  last  paragraph  clearly  indicate  that  the 
particular  state  from  which  an  unsteady  run  was  initiated  would  affect  the  subse¬ 
quent  flow  history.  In  fact,  more  generally  it  means  that  any  particular  unsteady 
state  is  dependent  on  the  previous  flow  history,  as  seemed  to  be  indicated  by  the 
deceleration  tests  with  the  uniform  tube. 
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SUMARY 

In  summary,  it  is  concluded  that  the  imposition  of  an  unsteady  transient 
produces  different  effects  for  the  two  basic  types  of  flow  investigated.,  as  follows* 

1.  For  cases  of  surface  resistance  caused  by  boundary  shear  stresses 

a)  With  acceleration  the  resistance  is  slightly  but  not  appreciably  greater  than 
for  the  equivalent  steady  state. 

b)  With  deceleration  the  resistance  is  appreciably  less  than  for  the  equivalent 
steady  state, 

c)  With  either  acceleration  or  deceleration,  it  appears  that  the  internal  flow 
structure  is  not  markedly  different  from  that  for  steady  states. 

2.  For  cases  of  form  type  resistance  associated  with  the  high  shear  and  generation 
and  diffusion  of  turbulence  accompanying  jet  formation 

a)  With  acceleration  the  resistance  is  appreciably  less  than  for  the  equivalent 
steady  state. 

b)  With  deceleration  the  resistance  is  appreciably  more  than  for  the  equivalent 
steady  state. 

c)  For  intense  jet  action  as  obtained  with  small  orifice  to  tube  diameter  ratios, 
it  appears  that  unsteadiness  produces  an  internal  flow  structure  that  is  no 
longer  comparable  to  any  steady  state  condition. 
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FIGURE  7  SMOOTH  CONDUIT  (AREA  RATIO  I.O) 


